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We calculate the ground-state energy of one and two-dimensional spatially inhomogeneous anti- 
ferromagnetic Heisenberg models for spins 1/2, 1, 3/2 and 2. Our calculations become possible as 
a consequence of the recent formulation of density-functional theory for Ifeisenberg models. The 
method is similar to spin-density-functional theory, but employs a local-density-type approximation 
designed specifically for the Heisenberg model, allowing us to explore parameter regimes that are 
hard to access by traditional methods, and to consider complications that are important specifically 
for nanomagnetic devices, such as the effects of impurities, finite-size, and boundary geometry, in 
chains, ladders, and higher-dimensional systems. 
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The study of low-dimensional spin systems is one of 
the central issues in the physics of correlated electrons. 
Model Hamiltonians of the Heisenberg type are widely 
used to study, e.g., antiferromagnetically coupled spin 
chains, ladders, and layers, ^'^ and constitute most use- 
ful models for strong correlations in undoped cuprates 
and manganites. Finite-size Heisenberg models, more- 
over, serve as paradigmatic models for the emerging field 
of nanomagnetism and spintronics. However, progress in 
the analysis of nanomagnetic devices requires capability 
to deal with 'real-life' complications, such as impurities of 
arbitrary size and location, boundaries of various shapes, 
and crossovers between finite and extended systems, or 
between one and higher dimensions. In the present paper 
we present a convenient and efficient numerical approach 
for calculating ground-state energies of antiferromagnetic 
Heisenberg models subjected to such complications. 

Traditional numerical methods are well suited for 
studying generic properties of homogeneous Heisenberg 
models, and have provided many important insights into 
the physics of magnetic systems. However, the very sig- 
nificant expenditure of computational effort required by 
methods such as group-theory-aided exact diagonaliza- 
tion. Quantum Monte Carlo (QMC),^ or the density- 
matrix renormalization group (DMRG),^ imposes rather 
strict limits on the size and complexity of systems that 
can be studied. The mean-field approximation can, in 
principle, be applied to systems of almost arbitrary size 
and complexity, but its neglect of correlation makes it un- 
reliable as a tool for studying many issues of current phys- 
ical interest in strongly correlated systems. In ah initio 
calculations, density-functional theory (DFT) provides a 
convenient and reliable way to go beyond the mean-field 
(Hartree) approximation, allowing to study real systems 
of considerable size and complexity.^'®'^ While DFT is 
in principle an exact reformulation of the many-body 
problem,^ its practical implementation requires the use 
of approximations for the exchange-correlation energy. 
Among the most important such approximations is the 
local-density approximation (LDA), the essence of which 



is to use the exchange-correlation energy of the uniform 
electron liquid locally as an approximation to the one of 
the crystal or molecule of interest.^'^ The formal frame- 
work of DFT and the LDA is not limited to the ah initio 
case, but can be applied to a large class of model Hamil- 
tonians, too.^'^°'"'^^'" 

In the present paper we study the ground-state en- 
ergy of one and two-dimensional antiferromagnetic sys- 
tems with broken translational symmetry, as a function 
of spin S and system size TV. Our choice of symmetry- 
breaking terms is dictated by complications expected to 
arise in real nanomagnetic devices, and includes bound- 
aries of diverse geometries, and different types of impuri- 
ties. Our calculations are performed within DFT for the 
Heisenberg model, 

H = J^S, -Sj+^B, -S,, (1) 

ij i 

where is a symmetry-breaking magnetic field, cou- 
pled to the spins. The simplest local-spin approxima- 
tion (LSA) for the correlation energy of Hamiltonian (1) 
reads^^ 

^i;lii[S,] = -0.36338 J^|S,| (2) 

i 

in one dimension, and 

£^Mf2[S.] = -0.316 J^IS.I (3) 

i 

in two dimensions. More sophisticated approximations 
are discussed in Ref. 12. Here |Si| is the modulus of the 
classical vector Si, the ground-state expectation value of 
the local spin operator S^. The total-energy functional 
then takes the form 

S[S,] ^JY,^.- S, + E^'^^iS,] + ^ • S„ (4) 

i-.j i 

where the first term on the right-hand side is the mean- 
field result, the second the LSA correlation correction to 
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it, and the last describes the couphng to (possibly spa- 
tially varying) magnetic fields. Below we set B,; = 
because we will mostly be interested in intrinsic inho- 
mogeneities such as surfaces and impurities, but both 
LSA functionals can be applied for nonzero = as 
well. Minimization of E[Si] with respect to Si yields the 
ground-state energy Eq. Minimization of the first (mean- 
field) term on its own results in the familiar antiferromag- 
nctic Nccl state with alternating up and down spins. For 
open boundary conditions^'* in one dimension this state 
has energy E^l^^{S) = -JS'^{N - 1), while for a two 
dimensional square lattice E^[^^{S) = -2JS'^N{N -1), 
where is the number of sites along the chain and along 
one side of the square, respectively. E^^-^ provides an 
additive correction to E^^ . 

First, we consider finite size Heisenberg chains. In 
Fig. 1 we compare, separately for each spin, the ground- 
state energies obtained in the mean-field approximation, 
in LSA, and by exact diagonalization. To be able to 
compare with exact results we have limited us to fairly 
small systems, but clearly mean-field and LSA calcula- 
tions can be performed easily for much larger systems. 
The LSA is seen to provide substantial improvement on 
the mean-field results, at negligible extra computational 
cost. However, we also observe that the LSA does not 
reproduce the oscillatory structure in the exact rc;siilts 
visible for small A^, but instead smoothly interpolates 
through it. As expected on physical grounds, the LSA 
becomes better as the system size increases. 

In these calculations the system is inhomogeneous only 
due to its finite size. However, the LSA concept can also 
be applied to spin Hamiltonians in which translational 
symmetry is broken more radically. Specifically, we now 
consider an impurity model of much current physical in- 
terest: an impurity of spin Si in a magnetic host ma- 
terial of spin S. This type of model is relevant, e.g., 
for the analysis of impurities in the antifcrromagnetic 
parent compounds of cuprate superconductors,^ impu- 
rities in correlated spin chains, and recent proposals 
for quantum computing.*®'^'' Fig. 2 shows the ground- 
state energy of a spin 1/2 antiferromagnetic Heisenberg 
chain, with open boundary conditions, in which one of 
the boundary spins has been replaced by a spin 1 ion, 
as schematically indicated in the inset. Mean-field, LSA, 
and exact values are given. The exact data show that 
for this type of impurity system the improvement on the 
mean-field approximation provided by the LSA is of the 
same order as in the homogeneous case studied in Fig. 1. 
For comparison purposes we have also included the LSA 
predictions for an impurity in the bulk, and for a sys- 
tem with two impurities, one in the bulk and one at the 
surface. Clearly, an antiferromagnetically coupled bulk 
impurity with Sj > S is energetically more favorable 
than a surface one, suggesting that self-assembled nano- 
magnetic systems with this type of impurity will tend 
to accumulate the impurities in the bulk, and leave the 
surface homogeneous. The ratio of the energy difference 
between bulk and surface impurities to temperature is 



the controlling parameter for studies of impurity migra- 
tion. We find that this parameter depends strongly on 
the impurity spin Si and system geometry. The extra 
computational effort for applying the LSA is small, and 
does not increase significantly if more than one impurity 
is present, thus opening the possibility to model realis- 
tic nanomagnetic systems, with hundreds or thousands 
of sites, with LSA. 

Next, we consider an antiferromagnetic background of 
spin 8=1/2 and varying size, with an impurity of spin 
Si > S at the surface. Fig. 3 shows how the combined 
system background+impurity approaches the thermody- 
namic limit for each value of spin Si. The lowest curve, 
representing the uniform system, approaches the ther- 
modynamic limit from below, and converges to it rather 
rapidly. The upper three curves represent the behaviour 
for different sizes of the impurity spin. The impurity 
curves approach the thermodynamic limit from above, at 
a rate that — due to the spin-dependence of the prefac- 
tors of the overall decay - decreases with increasing 
impurity spin. Physically, this behaviour can be under- 
stood by noting that for small N the impurity greatly 
increases the energy of the model as compared to the uni- 
form case, while, as the system size approaches infinity, 
the influence of the impurity localized at the boundary 
becomes negligible, and all curves converge to the value 
Eo/JN = -0.432, close to the value Eq/JN = -0.443 
obtained from the Bethe Ansatz for the uniform system. 

Next, we tiirn to the two-dimensional case. Here exact 
diagonalization becomes prohibitively expensive at even 
smaller system sizes than in one dimension, and alterna- 
tive methods, such as DMRG, also encounter very signif- 
icant difficulties. In Fig. 4 we plot the LSA ground-state 
energy of finite-size two-dimensional Heisenberg models. 
These are calculations without impurities, so that the in- 
homogeneity arises only from the presence of the system 
boundary, which breaks translational invariance. (As in 
one dimension, impurities of arbitrary size and location 
can be added to these systems without a significant in- 
crease in computational cost.) Given the LSA fmictional 
(3), the generation of the LSA data is a very simple ex- 
tension of the mean-field calculation. In spite of its sim- 
plicity, however, this calculation illustrates an important 
aspect of LSA: it has the same margin of error in higher 
dimensions as in one. This is rather untypical of many- 
body methods, which usually work better in certain di- 
mensionalities (often d = 1) than in others. 

Finally, we investigate an effect that does not exist in 
one dimension, namely the dependence of the ground- 
state energy on the geometry of the system. In Table I 
we compare the ground-state energies of five Heisenberg 
models with 100 sites, of shape 1 x 100, 2 x 50, 5 x 20, 
10 X 10 and 2 x 5 x 10. These systems represent, re- 
spectively, a one-dimensional chain of the type studied 
in Fig. 1, a spin ladder with constant coupling along legs 
and rungs, a multi-legged ladder, a square lattice of the 
type investigated in Fig. 4, and a representative three- 
dimensional structure. The data show that: (i) for larger 
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spins the effect of the system boundary becomes more 
pronounced. This may be relevant for the study of quan- 
tum corrals/'^ and differently shaped magnetic quantum 
dots, (ii) The higher-dimensional the system, the slower 
is the convergence to the thermodynamic limit: our result 
for the 100 site spin 1/2 chain differs only by 3% from the 
thermodynamical limit obtained from the Bethe-Ansatz, 
whereas with the same number of spins the energy of 
the 10 X 10 square lattice is about 9% from the thermo- 
dynamical limit reported in Ref. 20; a deviation of 3% 
from this limit is obtained only for a 50x50 lattice, (iii) 
The step from a chain to a two-legged ladder is much 
bigger than that from a two-legged ladder to a square 
lattice, suggesting that in terms of their ground-state 
energy even small ladders are more similar to extended 
two-dimensional systems than to chains, (iv) Compari- 
son of the columns 1 x 100 and 1 x 100(2(i) shows that 
simple extrapolation of a higher- dimensional functional 
to lower dimensionalities results in less negative energies 
than the correct functional. This observation may also 
be relevant for ab initio calculations applying LDA to 
low-dimensional systems. 

In summary, our results show that the Heisenberg LSA 
is quantitatively reliable, goes significantly beyond the 



mean-field approximation, and allows to calculate the 
energy of systems of considerable size and complexity. 
Similar results can be obtained with QMC,'^ but at orders 
of magnitude higher computational cost. (Unlike QMC, 
DFT is also not limited by a minus-sign problem.) QMC, 
DMRG and exact data provide benchmarks for testing 
the quality of improved density functional. On the other 
hand, LSA data on the effects of the boundary geome- 
try, the energetic influence of impurities of different sizes 
and locations, the approach to the thermodynamic limit, 
and effects of increased dimensionality may be useful in- 
formation for analysing real nanoscale antiferromagnetic 
systems, and can also serve as a guide for the application 
of alternative many-body methods to the same type of 
system. 
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FIG. 1: Ground-state energy of the finite amtiferromagnetic 
Heisenberg chain for various values of the spin S. For each 
spin, the diamonds represent the mean-field results, the cir- 
cles arc numerically precise benchmark values obtained in in- 
dependent many-body calculations,^** the full line is obtained 
with LSA, and the dashed horizontal line represents the value 
obtained in the thermodynamic limit, taten from Ref. 4. The 
LSA is seem to provide significant improvement on the mean- 
field values, for all S and N. 

FIG. 2: Ground-state energy of an antiferromagnetic S = 1/2 
Heisenberg chain with an ft = 1 impurity at the surface. 
Horizontal line: mean-field result. Dashed curve: LSA result. 
Circles: exact values, obtained by numerical diagonalization. 
Full curve: LSA result for same system, but with the impurity 
in the bulk. Dash-dotted curve: LSA results for same system, 
but with two impurities. 



FIG. 3: Ground-state energy of am antiferromagnetic S = 1/2 
Heisenberg chain with different impurities Si at the bound- 
ary. Mean-field values are not included because from Fig. 2 

we see that they are much inferior to LSA ones for this sys- 
tem. Crosses: homogeneous system. Diamonds: impurity 
spin Si — 1, open circles: impurity spin Si = 3/2, full cir- 
cles: impurity spin Si = 2. The arrow indicates the exact 
Bethe-Ansatz result Eq/NJ = -0.443147. 



FIG. 4: Ground-state energy of finite antiferromagnetic 
Heisenberg square lattices, as a function of spin and system 
size. Here N stands for the number of sites along the side 

of the square, hence the total number of sites is A'^'^. Dashed 
curves: mean-field (MF) results. Full curves: LSA results. 
Full circles: benchmark data obtained by exact diagonaliza- 
tion of small clusters. Arrow: N —> oo limit, from Ref. 20. 
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TABLE I: Ground-state energy —So/100 J of antiferromag- 
netic Heisenberg models of fixed size N = 100, and varying 
geometries and dimensionalities. For each dimensionality we 

used the proper functional, except in the column labeled 
1 X 100(2d), which contains values obtained for the Id system 
with the 2d functional (see discussion in main text). 

S 1 X 100 1 X 11M)(2(1) 2 X ")() 5 X 20 10 x 10 2 x ") x 10 

1/2 0.429 0.406 0.528 0.596 0.608 0.696 

1 1.353 1.306 1.796 2.066 2.116 2.491 
3/2 2.795 2.702 3.804 4.412 4.524 5.387 

2 4.687 4.592 6.552 7.632 7.832 9.382 
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